Pinholes may Mimic Tunneling 
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Interest in magnetic-tunnel junctions has prompted a re-examination of tunneling mea- 
surements through thin insulating films. In any study of metal- insulator- metal trilayers, one 
tries to eliminate the possibility of pinholes (small areas over which the thickness of the in- 
sulator goes to zero so that the upper and lower metals of the trilayer make direct contact). 
Recently, we have presented experimental evidence that ferromagnet-insulator-normal trilay- 
ers that appear from current- voltage plots to be pinhole-free may nonetheless in some cases 
harbor pinholes. Here, we show how pinholes may arise in a simple but realistic model of 
film deposition and that purely classical conduction through pinholes may mimic one aspect 
of tunneling, the exponential decay in current with insulating thickness. 
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Tunnel Junctions 



The construction of magnetic tunnel junctions with large room-temperature magnetore- 
sistance-*^ has triggered an intense research effort, as groups have applied them as magnetic- 
field sensors,^ memory devices, "^"^ and magnetic-medium read heads. ^ The current tech- 
nological drive toward a lower product of resistance and area,^~^ which implies the use of 
thinner and thinner insulating barriers, has re-opened the question of how to rule out the 
presence of pinholes, direct metal-metal contacts through the nominally insulating barrier. 
Recent high magnetoresistances of up to 300% at room temperature observed in magnetic 
nanocontacts^ raise the intriguing question of whether conduction through pinholes might 
actually contribute to the magnetoresistance of tunnel junctions. Similarly, any study of 
anomalous capacitance of magnetic tunnel junctions^"^^ must take pinholes into account. 
To exclude the possibility of pinholes, Rowell and others in the 1960s and 1970s developed 
a set of criteria to distinguish tunneling from other current paths. Four of these criteria 
continue to apply when neither electrode superconducts: (1) the exponential thickness de- 
pendence of the current, (2) the parabolic shape of the differential conductance as a function 
of voltage, (3) the scaling of the junction resistance with area, and (4) the temperature 
dependence of the conductivity. 

Wc have recently constructed a scries of junctions demonstrating experimentally that the 
second criterion applied alone cannot distinguish pinhole conduction from tunneling. ^^"^"^ 
The present work will show the first criterion similarly unreliable, since a purely classical 
conduction path through pinholes may under realistic assumptions mimic the exponential 
thickness dependence of the conduction current. This leaves only the temperature and area 
dependences of conductivity as good criteria for determining whether a junction contains 
pinholes. 

Kleinsasser et a/.^^ have shown that another of the proposed signatures of tunneling 
in superconducting-insulator- normal junctions, the subharmonic gap structure of Andreev 
reflection, may similarly reflect pinhole conduction, but the problem of flushing out pinholes 
from normal-metal junctions remains. 

Small pinholes may be invisible to surface microscopy, although recently "hot spots" 
have been observed using scanning tunneling microscopy. -"^^ Moreover, the large resistivity 
ratio between an insulator and a metal (about 10^-10^) requires that, over the junction area, 
pinhole regions of one part in 10^-10^ must be ruled out to ensure no pinhole conduction in 
parallel with tunneling. Generally speaking, this is very difficult to do and not commonly 
considered. Recognizing the technological, as well as scientific, importance of identifying 
pinholes in magnetic tunnel junctions, Schad et al. have developed a method for imaging 
pinholes through decoration by electrodeposited copper. This tool complements criterion 
(4) of temperature dependence and further highlights the risks of relying solely on thickness 
dependence and differential conductance. 

Tunneling versus Pinhole Conduction 

In a series of metal-insulator-mctal trilayers in which the insulating thickness, z, varies, 
the current at given voltage should decay exponentially in z (the applied voltage is assumed 
less than the insulating gap) . Such exponential decay has been cited as experimental evidence 
for good tunnel junctions. By the WKB approximation, the characteristic decay length 
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equals 



where m is an effective mass and a potential energy on the order of the barrier height. 

For 0=1/4 eV and m the bare electron mass, zq ~ 2A. Notably, for realistic parameters (f> 
and m, this decay length coincides roughly with the thickness of a single atomic layer. In 
contrast, a classical resistor supports a current inversely proportional to the thickness. 

The inverse relationship holds for a perfectly even layer, but real deposition processes 
leave an uneven insulator with possible pinholes. The simplest model for classical conduction 
by pinholes gives an exponential dependence of resistance on deposited thickness, mimicking 
quantum tunneling. Consider a metallic substrate on top of which we randomly deposit 
cubes of a perfect insulator to an average height /i (measured in monolayers). This is not 
a uniform height, so there may be pinholes. We then deposit a perfect conductor, making 
contact through any pinholes with the metal substrate. Since the insulator and overlayer are 
perfect, conduction is directly proportional to contact area and inversely to the metal- metal 
contact resistance, Rq. 

We consider deposition to take place on a regular two-dimensional grid, each cell of which 
may be occupied by any non-negative integral height (0, 1, . . .) of insulating particles. If 
each deposited particle can land randomly over any grid cell, the resulting heights follow a 
Poisson distribution. In the limit of an infinite two-dimensional grid, the probability that 
any given grid cell contains no insulating particles is exp(— //). For a large system, this is 
also the expected proportion of cells that will be unoccupied and so proportional to the 
conduction. Thus where tunneling can lead to an exponential-decay length of a monolayer 
for certain realistic parameters, the simple classical model always gives a decay length of one 
layer. 

This trivial model (perfect insulating blocks) predicts conductance that decays expo- 
nentially in coverage for all coverages. However, one expects a crossover to 1/// decay 
at larger coverages. The simplest extension of the model that might display such behavior 
deposits blocks randomly, but the insulators now have some finite resistance, R. In the 
regime of interest, R ^ Rq. As before, we deposit M blocks over an L x L lattice for a 
coverage /j, = M/ (with length measured in units of monolayer thickness). A face between 
the metal overlayer and an insulating block has resistance R, as does a face between two 
insulating blocks. Any connection to ground has contact resistance Rq. To make a simple 
model, we first turn off sideways conduction in the insulators, reducing the problem to that 
of an ensemble of independent columns of binomially-distributed resistors in parallel. We 
compare this "independent-column" model to numerical calculations below. Finally, we shall 
consider numerically an isotropic model that realistically models growth. 

Independent-column model 

Each independent column contains n resistors of value i?, where n is the column height, 
and one contact resistor of value Rq. Letting s — Rq/R, we have for the average column 
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conductance <(7>= \ {n+s)R / ' Poisson limit, — >■ oo, the scaled conductance 

oo 

S = R<a> = e-f Y . 

^ n\{n + s) 



n=0 --' - ■ (2) 

Jo * 
Here, iFi is a confluent hypergeometric function. 

One arrives at an asymptotic series for large /i through successive integrations by parts, 
exhibiting explicitly the approach to conductance: 



R<a> ~ ii ^ 



(3) 



where n is a cutoff that must be introduced for any finite /i. 

More interesting is the small-// limit, in which we recover the exponential decay of the 
trivial model, for multiplying the top equation in (2) by s, we have 

(00 n \ 

For small enough Rq <cr>~ e~^, representing exponential decay of conductance with 
coverage. 

We can estimate the crossover scale //q above which the decay ceases to be exponential by 
setting the "1" term in parentheses in (4) equal to the remainder. This yields the condition 

2 = iFi(s,s + l;^o) • (5) 
For s<0.5, we can replace iFi with the leading term in its asymptotic expansion, giving 

Ato « ln(2/s) . (6) 

The logarithmic form of (6) comes as something of a surprise: it means that the non-zero 
contact resistance cuts off the small-coverage, exponential regime no matter how tiny Rq is 
relative to the insulating R. 

The solid traces (for different s) in the semilog plot, figure 1, clearly show the exponential 
regime for small fi and the gradual deviation as pinhole conduction ceases to dominate. To 
display large- and small-// behaviors together, it is helpful to plot y = dlnS/d\nn = S'n/S 
against /i; the exponential regime is characterized by a constant slope, while in the asymp- 
totic l/ii regime, the graph approaches the constant y — —1 (inset to figure 1). In this 
independent- column model, the minimum of the graph corresponds closely to the estimate 
(6); later, in numerical simulations incorporating surface relaxation, we see deviations some- 
what earlier. 

We have assumed the ratio s small, as supported by the following estimate. Using 
Sharvin's semiclassical calculation-*^^ for contact resistance through a pinhole of diameter 
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FIGURE 1. Scaled conductances (solid traces) given by (2) for three values of the ratio s of contact to insulating 
resistance. According to the estimate (6), they begin to cross over from exponential to inverse-thickness behavior at 
about 5.3, 7.6, and 9.9. For comparison, wc define the heuristic measure /ii5 as the point at which the slope deviates 
from —1 by 15%; the corresponding estimates are 4.4, 7.2, and 9.8. The dotted lino (s = 0) is an exact exponential. 
The dot-daslicd line shows a numerical simulation (100 000 trials on 5 x 5 grids) for s = 0.0001 with J — —0.5, 
diffusion x = 1/3, and "up"-steps allowed. The inset plots the same curves in a way that makes explicit the approach 
to conductance (dlnS/dln/n = —1, horizontal line) as well as the initial exponential regime. Note that surface 
relaxation in this case causes a more rapid deviation from exponential decay (/U15 = 6.1) while also delaying the onset 
of conductance. 

~ lA, and assuming a typical metallic Fermi temperature and electronic density, we have 
contact resistance Rq ~ lO^^l. The "classical" resistance of an insulating block is less well 
defined, but a minimum resistivity ~ lO^Q cm suggests R>10^^Q, so that s — Ro/R< lO"-*^^, 
tiny indeed. We then estimate by (6) that insulating thicknesses up to about 24 monolayers 
should show exponential resistance. 

Numerical Simulations 

The preceding model provides a baseline for our simulations. In addition to restoring 
conduction through the sides of blocks, we adopt a model due to Pal and Landau^^"^^ to 
describe the motion of insulating blocks along the surface after they have fallen. Such models 
have found good agreement with experiment for long-wavelength and long-time (large-/x) 
features, although our interest lies in the opposite regime. 

For simphcity, we imagine the original metaUic layer to be flat.^^ After depositing some 

number of blocks in the manner of the previous section, wc allow a fraction of blocks at the 
interface to diffuse one unit along the surface. Wc employ a quasi-Mctropolis procedure to 
determine whether to accept or reject each move (the system is not at equilibrium) with 
an energy J for each face-to-face bond between blocks: caUing D the change in the number 
of bonds, the move is accepted unconditionally if it lowers the energy and otherwise with 
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probability exp{—DJ/kBT). Henceforth, we scale J by setting the thermal energy ksT — 1. 
Thus a negative J encourages clustering, while a positive J could represent stearic hindrance 
or an affinity for roughness. We adopt one more feature from the Pal-Landau models: blocks 
either may or may not jump "up" (further away from the substrate) while relaxing. 

We perform the simulation on an L x L grid, measuring conductance at set times, before 
averaging over many trials to reduce noise. Since long-wavclcngth features do not concern 
us in the present work, L can be as small as allowed by diffusion. Consider the diffusion 
parameter x defined as the fraction of interface blocks diffusing one lattice spacing per 
monolayer deposited. A simple random-walk argument establishes that for x of order unity, 
we can look at coverages as great as about 25 for a 5 x 5 grid. Numerically, we tested values 
of L from 5 to 60 and found no significant differences in conductances. The actual diffusion 
parameter x of course will depend on temperature and physisorption and chemisorption 
energies. We considered diffusion parameters x =0.33-10, roughly corresponding to a fast, 
but not unrealistic, deposition rate.^^ 

This model must behave the same as the independent columns in the two limits of zero 
and large coverage, but as figure 1 shows, surface relaxation can reduce the thickness at 
which deviations from exponential decay first become apparent, primarily by modifying the 
rate at which pinholes fill in. Because the deviation occurs continuously, we need a heuristic 
measure. The slope of the fine in a semilog plot (figure 1) determines the decay length; in 
both the simple model and the simulations, it always starts out at —1. Such a plot appears to 
differ by eye from a straight line when the slope changes by 15%; we call this crossover scale 
m^. As noted in the caption to figure 1, it also agrees reasonably well with the analytical 
crossover estimate (6) for the independent-column model with no relaxation. 

Figure 2 plots as a function of bond strength, J, for a few selected values of s and 

the relaxation parameter, x. We find that relaxation still permits a large value for fii^ for 
particular ranges of J. When blocks arc permitted to jump up, there is a tendency for 
pinholes to be extinguished less quickly than exponentially (see dot-dashed curve in figure 
1), since at moderate coverage more columns of height 1 turn into pinholes than vice-versa. 
A sufficiently attractive (negative) J counters this tendency by discouraging jumps out of 
height-1 columns, but too large a negative value of J leads to supcrcxponential decay. In the 
model without jumps, a sufficiently large positive (repulsive) J effectively cancels the effect 
of surface relaxation on pinholes. 

In some cases, we find that a reasonably large negative J can effectively alter the 
exponential-decay length over some range of coverages; see figure 3. To get a decay length 
substantially different from one monolayer, we need to modify the model. For example by 
prohibiting direct deposition on top of isolated pinholes (as by stearic hindrance), filling these 
only by surface relaxation, we compute both a longer decay length and a larger crossover 
scale (not shown). 

Implications for experiment 

We have concentrated on one signature of tunneling, showing that purely classical pinhole 
conduction can mimic the exponential dependence of resistance on barrier thickness. We 
have recently demonstrated experimentally that another signature, the nonlinear form of 
current I{V) (as a function of bias V")^^~^^ also may fail to distinguish classical conduction 
from tunneling. In a series of ferromagnet-insulator- metal junctions in which the "metal" 
is actually a superconductor, all samples could be fit well to a tunneling form above the 
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FIGURE 2. The horizontal axis measures the bond strength (negative attracts blocks), the vertical an estimated 
maximum thickness for the mimicry of tunneling, specifically the empirical measure fii^ at which the slope of a semilog 
plot, as in figure 1, deviates by 15% from —1. For comparison, the cutoff estimates (6) without relaxation for s = .01, 
.001, and .0001 arc 5.3, 7.6, and 9.9. Blocks may jump up in these traces. 

superconductor's transition temperature, Tc. However, some (but not all) showed Andreev 
reflection below T^, indicating the presence of pinhole shorts through the insulator. In 
other words, samples that appear (according to dl/dV and thickness dependence I{z)) to 
be good tunnel junctions may not be. 

We would expect that increasing temperature should enhance current through a tunnel 
junction (effectively lowering the barrier) while suppressing conduction through any metallic 
short. In fact, we found that the temperature dependence of current /(T) at zero bias does 
distinguish the pinhole-free from the shorted samples. 
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FIGURE 3. An example of a different effective slope: the dotted line is an exact exponential. Between the 
arrows (t ... j /i « 1...7), it is difficult to distinguish the simulated decay with relaxation (s = .0001, x = 1/3, 
J = —2) from a straight line of slope ~ —1.25, i.e., a decay length of 0.8 monolayers. No smoothing has been applied 
to the main graph. The inset shows the slope of the decay; the solid line applies parametric smoothing. 
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